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The EOS of strongly interacting matter at densities ten to fifteen orders of magnitude 
larger than the typical density of terrestrial macroscopic objects determines a number of 
neutron star properties, including the pattern of gravitational waves emitted following 
the excitation of nonradial oscillation modes. This paper reviews some of the approaches 
employed to model neutron star matter, as well as the prospects for obtaining new 
insights from the experimental study of gravitational waves emitted by neutron stars. 
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1. Introduction 

The equation of state (EOS) is a nontrivial relation linking the thermodynamic 
variables that specify the state of a physical systenffl The best known example is 
Boyle's ideal gas law, stating that the pressm-e of a collection of N noninteract- 
ing, pointlike classical particles, enclosed in a volume V , grows linearly with the 
temperature T and the average particle density n = N/V. 

The ideal gas law provides a good description of very dilute systems, in which 
interaction effects can be neglected. In general, the EOS can be written expanding 
the pressure, P, in powers of the density (throughout this paper, I will use units 
such that h = c = kB = 1, kB being Boltzmann's constant): 

P^nT [l + nB(T) + n^C{T) + ...] . (1) 

The coefficients appearing in the above series, called virial expansion, depend on 
temperature only, and describe the departure from the ideal gas law arising from 
interactions. The EOS carries a great deal of dynamical information and provides a 
link between measurable macroscopic quantities, such as pressure and temperature, 
and the forces acting between the constituents of the system at microscopic level. 

In this brief review, I will discuss theoretical models of strongly interacting mat- 
ter at densities up to fifteen orders of magnitude larger than the typical density 
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of terrestrial macroscopic objects. The EOS in this density regime largely deter- 
mines neutron star properties, including the pattern of gravitational waves emitted 
following the excitation of nonradial oscillation modes. 

Section 121 provides a summary of neutron star structure, while some of the ap- 
proaches employed to model neutron star matter are briefly reviewed in Section 
13 The dependence of the equilibrium properties of a nonrotating neutron star on 
the EOS describing matter in its interior is discussed in Sectional Section El is de- 
voted to the analysis of the imprint of the EOS on gravitational wave emission. The 
conclusions are stated in Section 

2. Overview of neutron star structure 

The internal structure of a neutron star, whose cross section is schematically illus- 
trated in Fig. ^ is believed to feature a sequence of layers of different composition. 




Fig. 1. Schematic illustration of a neutron star cross section. Note that the equilibrium density 
of uniform nuclear matter corresponds to ~ 2.7 X lO-'^^g/cm'^. 

The outer crust, about 300 m thick with density ranging from p ~ lO^g/cm'^ to 
the neutron drip density p^ — ^x lO^^g/cm , consists of a Coulomb lattice of heavy 
nuclei immersed in a degenerate electron gas. Moving from the surface toward the 
interior of the star the density increases, and so does the electron chemical potential. 
As a consequence, electron capture becomes more and more efficient and neutrons 
are produced in large number through weak interactions. 

At p = pd there are no more negative energy levels available to the neutrons, 
which are therefore forced to leak out of the nuclei. The inner crust, about 500 m 
thick and consisting of neutron rich nuclei immersed in a gas of electrons and neu- 
trons, sets in. Moving from its outer edge toward the center the density continues 
to increase, till nuclei start merging to give rise to structures of variable dimen- 
sionality, changing first from spheres into rods and eventually into slabs. Finally, at 
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/5 ~ 2 X lO^^g/cm'^, all structures disappear and neutron star matter reduces to a 
uniform fluid of neutrons protons and leptons in weak equilibrium. 

The density of the neutron star core ranges between ^ 2 x lO^"* g/cm'^, at 
the boundary with the inner crust, and a central value that can be as large as 
1 — 4 X 10^^ g/cm^. Just above 2 x IQ-'^'* g/cm'^ the ground state of matter is a 
uniform fluid of neutrons, protons and electrons. At any given density the fraction 
of protons, typically less than ~ 15%, is determined by the requirements of weak 
equilibrium and charge neutrality. At density slightly larger than the equilibrium 
density of uniform nuclear matter, po ^ 2.7 x lO^'^g/cm , the electron chemical 
potential exceeds the rest mass of the fi meson and the appearance of muons through 
the the process n — > p + /i~ + 17^ becomes energetically favorable. 

All models of EOS based on hadronic degrees of freedom predict that in the 
density range po < p < 2po neutron star matter consists mainly of neutrons, with 
the admixture of a small number of protons, electrons and muons. 

This picture may change significantly at larger density, with the appearance of 
strange baryons produced in weak interaction processes. For example, although the 
mass of the I]~ exceeds the neutron mass by more than 250 MeV, its production be- 
comes energetically allowed as soon as the sum of the neutron and electron chemical 
potentials equals the S~ chemical potential. 

Finally, as nucleons are known to be composite objects of size ^ 0.5 — 1.0 fm, 
corresponding to a density ~ lO^^g/cm , it is expected that, if the density in the 
neutron star core reaches this value, matter may undergo a transition to a new 
phase, in which quarks are no longer clustered into nucleons or hadrons. 

3. Models of neutron star matter EOS 

It has long been recognized that the Fermi gas model, leading to a simple poly- 
tropic EOS, predicts a maximum neutron star mass ^ 0.7 Mq, thus dramatically 
failing to explain the observed neutron star masse P. This failure clearly shows that 
neutron star equilibrium against gravitational colapse requires a pressure other than 
the degeneracy pressure, whose origin has to be traced back to hadronic interactions. 
Unfortunately, the need of including dynamical effects in the EOS is confronted with 
the complexity of the fundamental theory of strong interactions, quantum chromo 
dynamics (QCD). As a consequence, all available EOS of strongly interacting mat- 
ter have been obtained within models, based on the theoretical knowledge of the 
underlying dynamics and constrained, as much as possible, by experimental data. 

While the properties of matter in the outer crust can be obtained directly from 
nuclear data^, models of the EOS in the inner crust, corresponding to 4 x 10^^ < 
p < 2x lO^^'g/cm^, are somewhat based on extrapolations of the available empirical 
information, since the extremely neutron rich nuclei appearing in this density regime 
are not observed on eartlP. 

However, as most of the neutron star mass resides in the core region (typically, 
the crust olny accounts for less than ~ 2 % of the total mass), the star properties 
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discussed in this review are largely unaffected by the details of the EOS describing 
the crust. In what follows, I will only discuss the EOS of neutron star matter in the 
region of nuclear and supranuclear density, i.e. at p > po = 2.7 x lO^'^g/cm . 

Note that the results described in the following Sections have been obtained 
under the two standard assumptions, whose validity has long been establisheci^ 
that i) thermal effects can be disregarded and ii) neutron star matter is transparent 
to neutrinos produced in weak interaction processes. 

3.1. Nucleon matter in f3- equilibrium 

Models of the EOS of neutron star matter a,t po < p <2 — 3po rest on the premise 
that in this density regime nucleons can still be treated as individual particles. This 
assumption appears to be reasonable, as the shape of the proton charge distribution, 
obtained by Fourier transforming its measured elelctric form factor, is such that two 
protons separated by a distance of 1 fm hardly overlap one another. 

The EOS of nucleon matter is mainly obtained from two different approaches: 
nonrelativistic nuclear many-body theory (NMBT) and relativistic mean field theory 
(RMFT). 

In NMBT, nucleon matter is viewed as a collection of pointlike protons and 
neutrons, whose dynamics is described by the nonrelativistic Hamiltonian 

^ = E|i+E-^^+ E ^^'^^ (2) 

i j>i k>j>i 

where m and pi denote the nucleon mass and the momentum of the i-th particle, 
respectively, whereas Vij and Vijk describe two- and three-nucleon interactions. The 
two-nucleon potential, that reduces to the Yukawa one-pion-exchange potential at 
large distance, is obtained from a fit to the available data on the two-nucleon sys- 
tem, including both deuteron properties and ~ 4000 precisely determined nucleon- 
nucleon scattering phase shiftJ^. The purely phenomenological three-body term, 
Vijk Vij, is needed to reproduce the binding energies of the three-nucleon bound 
stateP and the empirical saturation properties of symmetric nuclear matte/^. 

The many-body Schrodinger equation associated with the hamiltonian of Eq. ^ 
can be solved exactly, using stochastic methods, for nuclei with mass number 
A < 10. The energies of the ground and low-lying excited states are in excellent 
agreement with the experimental datcPI. Accurate calculations can also be carried 
out for uniform nucleon matter, exploiting translational invariace and using eith er a 
variational approach based on cluster expansion and chain summation technique 
or G-matrix perturbation theorj^^. 

Within RMFT, based on the formalism of quantum field theory, nucleons are 
described as Dirac particles interacting through meson exchang In the simplest 
implement atio n of this approach the dynamics, modeled in terms of a scalar and a 
vector fielcP^, is described by the lagrangian density 

C = Cn + Cm + Ci7it , (3) 



I 



I 



5, 2008 6:57 WSPC/INSTRUCTION FILE benhar 



Neutron star matter equation of state and gravitational wave emission 5 

where 

Cn{x) = iPn{x) {i^ - m) ipNix) , (4) 
iPn{x) denotes the isospin doublet describing the proton and neutron fields and 

Cm{x) = -^F^,{x)F^^''ix) + ^mlV^{x)V'^ix) 

+ ^d^4>{x)d^c^{x) ~ \ml<j?{x) , (5) 

with Ff^u{x) — duVf^{x) ~ dfj_V^{x). The interaction term is defined in such a way as 
to give rise to a Yukawa-like meson exchange potential in the static limit: 

Cint{x) ^ ga(t>{x)'lpN{x)'lJjNix) - guVf^ix)'lpN{x)j'^1pN{x) ■ (6) 

In the above equations, Vf^{x) and <y{x) denote the vector and scalar meson fields, 
respectively, and m^^, rria, g^) and g^ are the corresponding masses and coupling 
constants. 

Unfortunately, the equations of motion obtained minimizing the action turn out 
to be tractable only in the mean field approximation, i.e. replacing the meson fields 
with their vacuum expectation values, which amounts to treating them as classical 
fields. Within this scheme the meson masses and coupling constants are determined 
by fitting the empirical properties of nuclear matter, i.e. binding energy, equilibrium 
density and compressibility. 

NMBT, while suffering from the obvious limitations inherent in its nonrelativis- 
tic nature, is strongly constrained by data and has been shown to possess a highly 
remarkable predictive power. On the other hand, RMFT is based on a very powerful 
and elegant formalism, but assumes a somewhat oversimplified dynamics, which is 
not constrained by nucleon-nucleon data. In addition, it is plagued by the uncer- 
tainty associated with the use of the mean field approximation, which is long known 
to fail in strongly correlated system^4i_ 

Using either NMBT or RMFT one can calculate the energy-density of nucleon 
matter at any baryon number density ns and proton fraction a;, as well as the 
proton and nucleon chemical potentials. 

The equilibrium conditions with respect to the processes (neutrinos are omitted, 
as they do not contribute to chemical equilibrium) 

n^p + e^,n^p + ^^, (7) 

obtained from the minimization of the total energy-density, e, subject the constraints 
of baryon number conservation and charge neutrality, are 

= + Me ; Me = M^i (8) 

and 

rip = Tie + ni_, , (9) 

where fii = {de/duijv is the chemical potential of the particle of type i {i = 
n,p, e, /i), whose number density is denoted by rii. Using the definition of the proton 
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fraction, x = Up/ns — ripjijin + rip), and exploiting charge neutrality, one can 
write all chemical potentials as a function of ub and x. Thus, for any given ns, x 
is uniquely determined by Eqs.ljHJ- 

Fig.|21shows the density dependence of the energy per baryon (left panel) and the 
proton, electron and muon fractions of /3-stable matter calculated by Akmal et al.-l^ 
using NMBT. For comparison, the left panel also includes the results corresponding 
to pure neutron matter (x = 0) and symmetric nuclear matter (x = 1/2). 




1.50 



Fig. 2. Left panel: energy per baryon of nucleon matter calculated by Akmal et alM^ as a function 
of baryon number density. The dashed, dot-dash and solid lines correspond to pure neutron matter, 
symmetric nuclear matter and /3-stable matter, respectively. The box represents the empirical 
equilibrium properties of symmetric nuclear matter. Right panel: density dep endence of the proton, 
electron and muon fractions of /3-stable matter obtained by Akmal et alM^ 



The EOS in the form P = P{e) can be readily obtained from the binding energy 
per baryon shown in Fig. [21 through 

and 

e = n,(A+^). (11) 
3.2. Strange hadronic matter 

As the density increases, different forms of matter, containing hadrons other than 
protons and neutrons, may become energetically favored. For example, the weak 
interaction process 

n + ^ T.^ + Ve , (12) 
leading to the appearance of a hyperon, sets in as soon as the condition 

Mn + A'e = Me- (13) 
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is fulfilled by the chemical potentials (typically at ub > 2no, no = .16 fm^'^ be- 
ing the baryon number density corresponding to equilibrium of symmetric nuclear 
matter). At larger density the production of A's is also energetically allowed. 

For any given ns, the relative aboundances of the different hadronic and lep- 
tonic species can be determined from the equations expressing the requirements of 
equilibrium with respect to weak interactions, conservation of baryon number and 
charge neutrality. 

In principle, both NMBT and RMFT can be generalized to take into account 
the appearance of hyperons. However, very little is k nown of their interactions. The 
available models of the hyperon-nucleon potential^ are only loosely constrained by 
few data, while no empirical information is available on hyperon-hyperon interac- 
tions. 




Fig. 3. Composition of neutron star matter calculated by Baldo et allilil using G-matrix pertur- 
bation theory. Top panel: nucleons and leptons only. Middle panel: nucleons, leptons and noninter- 
acting hyperons (S~ and A). Lower panel: hyperon-nucleon interaction included. The difference 
between the solid and dashed lines provides a measure of the model dependence associated with 
the NN potential Vij (see Eq. j^). 



Fig. 13 shows the composition of neutron star matter resulting from the NMBT 
calculations of Baldo et Comparison between the middle and lower panel 

illustrates the r ole of hyperon-nucleon interactions, modeled using the potential 
of Maessen et alU^, whose inclusion pushes the threshold of hyperon production 
towards larger density. Hyperon-hyperon interactions are neglected altogether. 

The transition from nucleon matter to hadronic matter makes the EOS softer, 
thus leading to a larger compressibility. This feature can be easily understoood, as 
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processes like the one of Eq. (|12|) replace particles carrying large Fermi energies 
with more dilute, and therefore less energetic, strange baryons. 

3.3. Quark matter 

Due to the complexity of QCD, a first principle description of the EOS of quark 
matter at high density and zero temperature is out of reach of the ex istin g compu- 
tational approaches. Following the pioneering work of Baym and Chir^l a number 
of authors have car ried out numerical studies of quark matter based on the simple 
MIT bag modePI. Within this model the main features of QCD, namely confine- 
ment and asymptotic freedom, are implemented through the assumptions that: i) 
quarks occur in color neutral clusters confined to a finite region of space (the bag), 
whose volume is limited by the pressure of the QCD vacuum (the bag constant B), 
and ii) residual interactions between quarks are weak, and can be treated in low 
order perturbation theory. 

Neglecting quark masses, the bag model EOS, at first order in the color coupling 
constant a,,, can be obtained in closed form from the relations linking pressure and 
energy-density to the quark chemical potentials 

--i(i-^)i:.j--. (") 

and 

^ ^ / 

where the sum runs over the active flavors, typically u, d and s, and fJ-f is the 
chemical potential of the quark of flavor /. 

For any baryon desity, quark densities are dictated by the requirements of baryon 
number conservation, charge neutrality and weak equilibrium. In the case of three 
active flavors one finds 

1 2 11 

nB ^ -^{riu + nd + ns) , -Uu - -Ud - -Us - Ue = (16) 

and 

^J-d = fJ-s = fJ-u + Me ■ (17) 

Note that in the above equations the possible appearance of /x mesons is not taken 
into account, as the results of numerical calculations show that the electron chemical 
potential never exceeds the muon mass. 

Fig.^shows the composition of charge neutral quark matter in weak equilibrium, 
obtained from the MIT bag model. The calculations have been carried out including 
massless u and d quarks and strange quarks of mass rus = 150 MeV, and setting 
B = 200 MeV/fm"^ and ~ 0.5. It appears that at large densities quarks of the 
three different flavors are present in equal number, and leptons are no longer needed 
to guarantee charge neutrality. 
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Fig. 4. Composition of charge neutral matter of u, d and s quarks and electrons in weak equi- 
librium, obtained from the MIT bag model setting mu = = 0, rris = 150 MeV, B = 200 
MeV/fm^ and cxs = 0.5. 



4. EOS and properties of nonrotating neutron stars 

Plugging 1 
o 

equations^ 



Plugging the EOS, P = P{e), into the Tolman Oppenheimer Volkoff (TOV) 
.202 



dPjr) _ [e{r) + P(r)] [Af (r) + ATrr^Pjr)] 
dr r2[l-2GM(r)/r] ' ^ ' 

where G denotes the gravitational constant, and 

M{r) = 47r / r'^dre{r) , (19) 



one can obtain the properties of stable nonrotating neutron stars. Eqs. H18|l and 
(|19|l are solved by integrating outwards with the initial condition e(r = 0) = Ec- For 
any given value of the central density, Ec, the star radius R is determined by the 
condition P{R) = and its mass M = M{R) is given by Eq. fH^ . 

Plotted as a function of central density, the neutron star mass exhibits a maxi- 
mum, whose value M^ax is mostly determined by the stiffness of the EOS. Stiffer 
EOS, corresponding to more incompressible neutron star matter, lead to larger 
J^max- Therefore, in principle, measurements of neutron star masses may be used 
to constrain the models of EOS at p > po. 

Unfortunately, comparison of the calculated Mmax with the neutron star masses 
obtained from observations, ranging between ~ 1.1 and '^1.9 Mgpl^, does not pro- 
vide a stringent test on the EOS, as most models predict a stable star configuration 
with mass compatible with the data. A stronger constraint may soon come fr om 
measurements of the neutron star mass-radius ratio. It has been recently reported^ 
that the Iron and Oxygen transitions observed in the spectra of 28 bursts of the 
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X-ray binary EXO0748-676 correspond to a gravitational redshift z — 0.35, yielding 
in turn a mass-radius ratio of the source M/R — 0.153 A/Q/km. 

The results of Cottam et al.^2 are still somewhat controversial and need to be 
confirmed. However, Fig. [3 shows that the M{R) relations corresponding to EOS 
obtained from NMBT including only nucleon degrees of freedom are consistent with 
the redshift measurement. The possible appearance of deconfined quark matte r in 
a small region in the center of the star does not dramatically change the picturJ^, 
while the occurrence of a transition to hyperonic matter at densities as low as twice 
the equilibrium density of nuclear matter seems to be ruled out^^. From Fig.^it also 
appears that none of the considered EOS, all of them obtained from nonrelativistic 
NMBT, predicts a M{R) curve extending into the region, forbidden by causality, in 

no 

which the speed of sound in matter exceeds the speed of light^-. 




Fig. 5. Mass r adiu s relation for diflfercnt EOS. Solid lin e: nu cleons onl}lii4 Dotdash: nucle- 
ons and liyperonjl^ Dashes: nucleons and deconfined quarks,-'^!. The horizontal lines denote the 
observational limits on the neutron star mass, wh erea s the other straight lines correspond to the 
gravitational redshift measurement of Cottam et all^^and to the boundary of the region forbidden 
by causality. 



5. Gravitational waves from neutron stars 

When a neutron star is perturbed by some external or internal event, it can be 
set into non radial oscillations, emitting gravitational waves at the charateristic 
frequencies of its quasi-normal modes. This may happen, for example, as a conse- 
quence of a glitch, a close interaction with an orbital companion, a phase transition 
occurring in its inner core or in the aftermath of a gravitational collapse. The fre- 
quencies and damping times of the quasi-normal modes carry information on the 
structure of the star and the properties of matter in its interior. 

Quasi-normal modes are classified according to the source of the restoring force 
which prevails in bringing the perturbed element of fluid back to the equilibrium po- 
sition. Thus, we have a g-mode if the restoring force is mainly provided by buoyancy, 
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or a p-mode if it is due to a gradient of pressure. The frequencies of the g-modes 
are lower than those of the p-modes, the two sets being separated by the frequency 
of the fundamental f-mode, associated with global oscillations of the fluid. General 
relativity also predicts the existence of additional mode s, calle d w-modes, which are 
purely gravitational, as they do not induce fluid motio n^ l ^ l The w-modes can be 
both polar (even parity) and axial (odd parity). They are highly damped and, in 
general, their frequencies are higher than the p-mode frequencies. 

To see an example of how the oscillation frequencies of quasi-normal modes 
depend on the neutron star matter EOS, consider the axial w-modes of a nonrotating 
star. In this case, the complex frequencies are eigenvalues of a Schroding er-li ke 
equation whose potential Ve{r) explicitely depend upon the EOS according tcP^ 

Vi{r) = —3- {i{e + l)r + [e(r) - P(r)] - 6M{r)} , (20) 

where 

diy 1 dP 



(21) 

dr [e(r) + P{r)] dr ^ ' 

Using a set of different EOS one obtains strongly damped eigenmoded^, whose 
frequencies exhibit the behavior displayed in Fig. |3 
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Fig. 6. Frequencies of the first axial w-mode of a nonrotating neutron star, plotte d as a function 
of the star compactness. The different curves correspond to different models of EOsCHj 



The pattern emerging from Fig. El strictly reflects the stiffness of the different 
EOS in the relevant density region (typically po < p < 5po), softer EOS corre- 
ponding to higher frequencies. For example, the curve labelled B h as been obtained 
from a model based on NMBT and including strange baryonJ22l which leads to a 
very soft EOS and a maximum mass of ~ 1.4 Mq. On the other hand, the curve 
labelled L corresponds to a RMFT calculation including nucleons onlj*2fli yielding 
a maximum mass of ~ 2.7 Mq. 

In 1998, Andersson and KokkotaJ^ computed the frequencies of the f-mode, 
the first p-mode and the first polar w-mode of a non rotating neutron star for a 
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number of EOS available at that time. They fitted the results of their calculations 
with appropriate universal functions of the macroscopic properties of the star, i.e. 
mass and radius, and showed how these empirical relations could be used to put 
constraints on these quantities if the frequency of one or more modes could be 
identified in a detected gravitational signal. 

A similar analysis has been recently carried out by Benhar et alP^, who have 
included in their fit results obtained from state-of-the art EOS. As an example. 
Fig. 13 shows the frequencies of the fundamental mode corresponding to different 
EOS, plotted clS 3jS 3b function of the square root of the average density. It appears 
that the frequencies corresponding to different EOS tend to scale to a straight line. 
The sizable displacement (~ 100 Hz) of the scaling line of Benhar et al. (thick 
solid straight line), with respect to the one obtained by Andersson and KokkotaJ^ 
(dashed straight line), reflects more than a decade of improvements of theoretical 
models of the EOS. 

The damping time of the f-mode and the frequencies and damping times of the 
p- and w-modes also exhibit scaling, when plotted as a function of th e comp actness 
(M/R). As a consequence, as pointed out by Anderson and Kokkota d'^ l ^"^ !, identi- 
fication of the frequency and damping time of some of these modes in a detected 
gravitational signal would allow one to determine the radius of the source knowing 
its mass. 



f-mode frequency 



X 




Fig. 7 

for different E09^^ The fits of Andersson and Kol^kota; 
fit and New fit, respectively. 



Frequency of tlie f-mode plotted as a function of th e square root of the average density 

i^and Benhar et ah^are labelled AK 



In a ddition to carrying out the analysis proposed by Anderson and 
KokkotaJSIEni, one may want to address a different, more direct, question. Will 
the simultaneous knowledge of the mode frequencies and the mass of the star, its 
only accurately measured observable, provide a severe test to discriminate among 
different EOS? 

As an example, let us consider the fundamental mode. Numerical simulations 
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show that this is the mode which is mostly hkely to be excited in many astrophysical 
processes and its damping time is quite long, so that it should appear as a sharp 
peak in the spectrum of the gravitational signal. 

Figure |H1 shows the mass dependence of the f-mode freq uencies obtained from 
the EOS models included in the analysis of Benhar et alP^. 

Comparison between the curve labelled APR2 and those labelled APRB120 and 
APRB200 indicates that the presence of quark matter in the star inner core does 
not significantly affect the pulsation properties of the star. This is a general feature, 
also observed in the behavior of p- and w-modes. 

The BBSl and APR2 EOS, based on similar dynamical models, turn out to yield 
appreciably different f-mode frequencies. This is likely to be ascribed to the effect of 
the relativistic corrections included in the APR2 EOS and to different treatments 
of three-nucleon interactions. 

The transition to hyperonic matter, predicted by the BBS2 model, produces a 
sizable softening of the EOS, thus leading to stable neutron star configurations of 
very low mass. As a consequence, the corresponding f-mode frequency is significantly 
higher than those obtained from the other models. So much higher, in fact, that its 
detection would provide a clear signature of the presence of hyperons in the neutron 
star core. 



f-mode frequency 
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Fig. 8. Frequency of the f-mode, plotted as a function of the mass of the emitting star for different 
models of EO 
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It is also interesting to compare the f-mode frequencies corresponding to mod- 
els BBS2 and G240, as they both predict the occurrence of strange baryons but 
are obtained from different theoretical approaches, NMBT and RMFT, based on 
different descriptions of the underlying dynamics. The behavior oi displayed in 
Fig. |S| directly reflects the relations between mass and central density obtained from 
the two EOS, larger frequencies being always associated with larger densities. For 
example, the configurations of mass ^ 1.2 Mq obtained from the G240 and BSS2 
have central densities ~ 7x10^'' g/cm'^ and ^ 2.5x10^^ g/cm^, respectively. On the 
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other hand, the G240 model requires a central density of ^ 2.5x10^'^ g/cm^ to reach 
a mass of ^ 1.55 Mq and a consequent i/f equal to that of the BBS2 model. 

The curves labelled SSI and SS2 correspond to strange stars, entirely made 
of quark matter with equal number of up, down and strange quarks. The peculiar 
properties of these stars, which are also a pparen t in the mass-radius relation, largely 
depend on their being self-bound object; 



6. Conclusions 

The results briefly reviewed in Section^lsuggest that the observation of gravitational 
waves emitted by neutron stars may provide valuable new insight, allowing one to 
severely constrain theoretical models of strongly interacting matter at high density 
and zero temperature. 

To address the question of whether the existing gravitational antennas may 
actually be able to detect these signals let us consider, as an example, a neutron 
star of mass M = 1.4 Mq described by the EOS of Akmal et al.^^, and assume that 
its fundamental mode (whose frequency and damping time are Vf = 1983 Hz and 
Tf = 0.184 s, respectiveljE^) has been excited by some external or internal event. 
The signal emitted by the star can be modeled as a damped sinusoi 

h{t) = ^e^*"'-*)/^^ sin [2T:vf [t - tarr)] , (22) 

where iarr is the arrival time and A is the mode amplitude. The energy stored into 
the mode can be estimated by integrating the energy flux 

dEmodc 2 I \ |2 Coo^ 

where h^iy) is the Fourier transform of h(t), over surface (S) and frequency (v). 

The signal to noise ratio (SNR) can be expressed in terms of h{i') and the 
noise power spectral density of the detector, S'„(i^). In the case of the ground based 
interferometric antenna VIRGcE3 one finds that SNR=5 corresponds to _E/ ~ 
6 X 10^^ A/0 for a source in our Galaxy (distance from Earth ^ 10 kpc) and 
Ef ^ 1.3 A/0 for a source in the VIRGO cluster (d ~ 15 Mpc). Similar estimates 
can be found for the detectors LIGO, GEO and TAMA. 

These numbers suggest that it is unlikely that the first generation of interfero- 
metric antennas will detect gravitational waves emitted by an oscillating neutron 
star. However, detection will become possible with the next generation of detectors, 
which are expected to be much more sensitive at frequencies above 1-2 kHz. 

Identification of the modes in a detected signal, combined with the knowledge 
of the mass of the emitting star will provide critical new information on its internal 
structure, thus opening the era of gravitational wave asteroseismology. 
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The results discussed in Sections 01 and have been obtained in collaboration with 
E. Berti, V. Ferrari, L. Gualtieri and R. Rubino. 
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